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I.  INTRODUCTION 


Establishing  the  structure  and  bonding  of  borane  and  carborane  molecules 
has  been  a  significant  area  of  chemical  research  for  many  decades.  While  much 
has  been  learned,  as  Beaudet's  recent  review  emphasizes,  there  are  still  many 
borane  and  carborane  systems  for  which  the  molecular  structure  is  poorly 
resolved  or  unknown.  It  is  therefore  not  surprising  that  there  is  little 
reliable  empirical  data  describing  the  basic  thermochemical  properties  of  the 
boranes.  We  present  here  a  discussion  of  structures,  harmonic  vibrational 
frequencies,  and  heats  of  formation  for  a  series  of  small  boranes.  The 
predictions  of  molecular  heats  of  formation  presented  here  use  a  scheme  which 
is  similar  to  the  bond-additivity  correction  (BAC)  method  pioneered  by  Melius 
and  Binkley. ^  The  technique  is  based  on  the  observation  that  although  the 
accurate  first  principles  calculation  of  bond  dissociation  energies  is 
prohibitive  for  most  systems  of  interest,  the  error  incurred  using  a 
correlated  quantum  chemical  method  with  a  modest  basis  set  is  fairly  constant 
for  a  given  bond  type.  This  notion  has  been  eKploited  by  Melius  and  Binkley 
by  combining  quantum  chemical  calculations  and  bond-additive  correction 
factors  to  calculate  the  heats  of  formation  for  hundreds  of  molecules,  in  manv 
cases  providing  important  data  that  is  not  experimentally  accessible. 

The  success  of  the  RAC  technique  is  due  in  part  to  the  availability  of 
accurate  experimental  heats  of  formation  for  a  number  of  the  molecules  in 
question.  This  data  is  compared  with  the  calculated  results  to  determine 
correction  factors  for  each  type  of  chemical  bond. 

The  BAC  technique,  as  described,  cannot  be  applied  to  the  boranes  because 
there  are  not  nearly  enough  known  heats  of  formation  to  determine  the  bond 
correction  factors.  We  seek  to  circumvent  this  problem  by  using  the 
capabilities  of  modern  quantum  chemical  methods  to  predict  acceptably  accurate 
atomization  energies  for  certain  small  prototypical  molecules.  With  this 
information  and  available  atomic  heats  of  formation,  we  can  directly  compute 
accurate  values  of  heats  of  formation  for  these  species.  These  calculations 
are  at  a  high  level  and  cannot  be  applied  generally.  However,  by  comparing 
these  high  level  calculations  to  calculations  performed  on  the  same  molecules 
at  a  more  widely  applicable  level  of  theory,  we  determine  bond  correction 
factors . 


We  report  here  results  for  11  molecules  containing  one  or  two  boron  atoms 
and  between  zero  and  six  hydrogen  atoms.  Many  of  these  compounds  have  been 
studied  at  both  levels  of  theory  alluded  to  above.  Bond  correction  factors 
are  determined  for  terminal  BH  bonds,  three  center  BHB  bonds,  and  BB  bonds. 
Structures,  harmonic  vibrational  frequencies,  and  heats  of  formation  are 
presented  for  all  of  these  species. 

II.  DESCRIPTION  OF  THE  CALCULATIONS 


Molecular  structures  and  harmonic  frequencies  for  all  the  compounds 
studied  were  obtained  using  the  Hartree-Fock  method  with  a  double  zeta  plus 
(  DZP)  polarization  basis  set  The  Gaussi  an-B?.  computer  codes  were  used  for 

these  studies.  Standard  double-zeta  basis  sets  were  augmented  by 
polarization  functions  with  exponents  optimized  for  correlation  energy 
calculations^  (0.421  for  the  boron  and  0.72S  for  hydrogen).  The  SCF  optimized 
geometries  were  used  to  determine  atomization  energies  by  computing 
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correlation  energy  using  full  fourth-order 
in  Gaussian-82 .  Most  of  these  calculations 
polarization  basis.  Several  calculations, 
set  built  upon  the  6-11 1G  basis  set.'7  The 
functions  in  the  valence  shell,  three  sets 
polarization  space  and  a  set  of  functions  1 
space.  Exponents  for  these  functions  were 
al.  The  atomization  energy  is  computed  as 
atomic  electronic  energies  and  the  electron 


perturbation  theory  as  implemented 
employed  the  double  zeta  plus 
however,  employed  an  extended  basis 
extended  basis  added  diffuse 
of  functions  in  each  atom's  first 
n  each  atom's  second  polarization 
those  recommended  by  Frisch,  et 
the  difference  between  the  sume  of 
ic  energy  of  the  molecule. 


III.  RESULTS  FOR  MONOBORANES 


Pople  and  coworkers  have  reported  accurate 
three  B,  monoboranes  using  an  extended  basis  set 

i  Q 

perturbation  theory  calculations.  Our  calculat 
reported  in  several  ways.  Rather  than  including 
extensions  as  seperate  contributions  to  the  elec 
fourth-order  energy  in  a  single  calculation  that 
functions.  Since  we  used  the  full  basis  set,  we 
energy  contributions  from  higher-order  perturbat 
description  of  our  methods  is  given  in  an  earlie 
report  fourth-order  results  using  both  the  doubl 
set,  and  the  extended  basis  set. 


atomization  energies  for  the 
in  full  four th-order 
ions  differ  from  those 
the  effects  of  the  basis  set 
tronic  energy,  we  computed  the 
included  all  the  basis 
did  not  include  estimates  of 
ion  theory  terms.  A  thorough 
r  contribution.^  Here  we 
e-zeta  plus  polarization  basis 


Table  1  reports  the  structural  parameters  and  scaled  harmonic  frequencies 
for  BH,  Bl^  and  BH^  and  Table  2  summarizes  the  perturbation  theory  energy 
predictions  for  these  compounds  and  both  the  hydrogen  and  boron  atoms. 


Table  1.  Structural  Parameters  and  Frequencies:  BH,  BH^ ,  and  BH^ 


Harmonic 


kbh 

(A) 

^  HBH 

(  Deg  . ) 

Frequencies 

(cm  1 ) 

Molecule 

Theory 

Exp . 

Theory 

Exp . 

Theory 

Exp . 

BH 

1  .218 

1  .214 

— 

2474 

2368 

BH? 

I  .190 

1  .178* 

1  28  .0 

180  .0* 

2841 

2650 

2681 

1430 

1076 

840(2) 

BHo 

l  .195 

1  .164* 

120  .0 

1 20  .0* 

2776(2) 

2976(2) 

1 26  3(2)  1765(2) 


*Data  estimates  from  Reference  11. 
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Table  2.  Total  and  Atomization  Energies:  H,  B,  BH ,  BH2  ,  and  BH3 


Molecule 
( Atom) 

H 

MP4 

DZP 

-0  ,497b 

Energy  (h)* 
Extended 

-0  .4998 

Atomization 

DZP 

Energy  ( h) 
Extended 

B 

-24  .5993 

-24  .6117 

BH 

-25  .2243 

-25  .2448 

0  .1274 

0 .1336 

bh2 

-25  .8559 

-25  .8819 

0  .2614 

0  .2709 

BH-. 

-2b  .5259 

-26  .5578 

0  .4338 

0  .4467 

-2b  .53676( FC) 


*F,nergy  units  are  hartrees.  One  hartree  =  627  .5  kcal/mole. 


Of  all  the  b°r<'n  hydrides  discussed  in  this  paper,  the  BH  molecule 

appears  to  have  the  most  accurately  known  structure  and  heat  of  formation. 

r  1  4*  J  12 

The  X  E  state  has  been  studied  using  spectroscopic  techniques  by  Gaydon. 

There  is  excellent  agreement  between  the  experimental  and  theoretical  values 
for  the  equilibrium  bond  length,  while  the  vibrational  frequency  prediction  is 
too  high.  It  should  he  noted  that  the  bond  length  for  the  diatomic  is 
substantially  longer  than  the  bond  lengths  for  the  other  Rj  boranes.  The 
electronic  energy  calculations  display  a  large  energy  difference  between  the 
two  basis  sets,  more  than  20  mi l 1 ihar t rees  f  20  mi  1 1 ihart rees  corresponds  to 
ca.  12. 5  kcal/mole).  The  predicted  atomization  energy  differences,  on  the 
other  hand,  are  relatively  small,  ca.  3.6  kcal/mole.  Were  it  assumed  that  the 
extended  basis  set  calculation  yielded  the  correct  atomization  energy,  then 
the  calculated  bond-additivity  correction  relevant  to  the  DZP  basis  for  this 
molecule  is  equal  to  that  difference.  Since  empirically  derived  BACs  have 
been  found  to  depend  parametrically  on  bond  length,  this  value  provides  a  long 
bond  estimate  of  the  correction  factor.  Our  best  prediction  of  the 
atomization  energy  is  82.7  kcal/mole,  a  number  that  agrees  well  with  Pople’s 
prediction  of  85.3  kcal/mole  and  with  the  experimental ly  derived  value, 

82  .2  kcal/mole  . 

In  contrast  to  the  BH  molecule,  there  is  no  reliable  empirical  data  for 
the  RH2  molecule.  The  estimates  of  the  structure,  vibrational  frequencies  and 
heat  of  formation  presented  in  the  JANNAF  tables  are  incorrect.  The  predicted 
structural  parameters  presented  in  Table  2  are  those  for  the  lowest  energy 
doublet  state  of  this  free  radical  species.  The  spin-eigenfunc.tion  of  the  i'HF 
wavefunction  is  almost  equal  to  two,  so  we  expect  that  the  effects  of  "spin- 
contamination"  are  small.  The  total  energy  difference  between  the  DZP  and 
extended  basis  set  calculations  is  large,  but  as  was  the  case  with  BH ,  the 
predicted  atomization  energies  differ  by  a  relatively  small  amount,  .0095  h. 
This  is  equivalent  to  almost  exactly  6  kcal/mole  difference  between  the  two 
atomization  energy  predictions. 
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The  last  of  the  monoborane  compounds  is  BH^  .  While  one  expects  to  find  a 
significant  literature  detailing  the  molecular  structure  and  properties  of 
borane[3],  there  is  surprisingly  little  data  available  for  this  molecule  (the 
notation  borane[3)  denotes  a  single  boron  atom  bonded  to  these  hydrogen 
atoms).  There  are  no  reliable  measurements  of  the  structure  of  the  molecule 
nor  of  the  vibrational  frequencies.  The  published  heat  of  formation  values 
are  all  derived  indirectly  based  upon  the  value  assumed  for  diborane,  ^2^6* 

The  molecular  structure  parameters  and  vibrational  frequency  estimates 
reported  here  agree  well  with  the  values  predicted  in  other  theoretical 
studies  of  the  molecule.^ In  this  case,  the  electronic  energy  difference 
between  our  two  basis  sets  is  19  mil lihartrees ,  while  the  predicted 
atomization  energy  difference  is  almost  13  mi L 1 ihart rees,  or  8.1  kcal/mole. 


As  the  number  of  bonds  increases,  the  dl 


^tween  the  atomization 


energy  computed  with  the  extended  basis  set  differs  more  and  more  from  that 
computed  with  the  standard  basis,  while  the  total  electronic  energy  difference 
for  the  molecules  remains  relatively  stable.  This  dichotomy  exemplifies  two 
effects.  The  total  energy  difference  between  the  two  basis  sets  is 
principally  an  atomic  effect.  In  the  case  of  these  three  compounds,  the 
contribution  of  the  boron  atom  dominates.  On  the  other  hand,  the  increasing 
disparity  in  the  atomization  energy  results  with  the  increase  in  the  number  of 
R-H  bonds  clearly  demonstrates  the  importance  of  the  extended  basis  set  in 
obtaining  accurate  predictions  for  this  quantity.  As  the  molecular  structure 
becomes  more  complicated,  the  magnitude  of  the  error  that  must  be  corrected 
grows  larger.  That  the  difference  depends  upon  the  number  of  bonds  indicates 
that  a  bond-additivity  correction  scheme  may  provide  a  method  for 
systematically  estimating  thermochemical  parameters  for  the  boranes. 


Overall  these  data  compare  reasonably  well  with  Pople's  results,  in  spite 
of  significant  differences  in  our  application  of  the  perturbation  theory 
methods.  The  predicted  data  are  certainly  better  than  the  data  available  from 
experimental  sources.  This  is  encouraging,  since  there  is  virtually  no 
experimental  data  for  any  of  the  diborane  species. 


IV.  RESULTS  FOR  1)  I  BO  RANKS 


we  '..scribe  results  for  ni  ne  di  boranes  and  the  lowest  closed  shell 
singlet  state  of  the  B ^  molecule.  Except  for  the  diborane[6]  molecule,  we 
report  electronic  energy  results  only  for  the  double  zeta  plus  polarization 
basis  set.  Schematic  structures  for  the  molecules  are  presented  in  the 
figures.  Detailed  structural  parameters  are  summarized  in  the  tables.  Only 
boron-hydrogen  bonds  were  encountered  In  the  monoborane  calculations.  The 
diboranes  provide  data  on  both  boron-boron  sigma  and  pi  bonds,  as  well  as 
examples  of  the  characteristic  BHB  three-center  bond.  Oata  on  these  bonding 
tvpes  are  essential  to  the  concept  of  bond-related  correction  factors. 

Table  3  contains  data  on  the  B j  molecule.  We  are  interested  in  the 
closed  shell  state  of  this  molecule.  The  bond  length  is  predicted  to  he  long, 
almost  two  Angstroms,  but  the  second  derivative  calculation  identifies  the 
extremum  as  a  stable  point.  We  have  computed  the  atomization  energy  for  this 
structure  using  the  extended  basis  set,  and  the  molecule  is  predicted  to  be 
bound.  There  is  the  possibility  that  our  molecular  structure  prediction  is, 
in  this  case,  not  as  accurate  as  those  obtained  in  the  case  of  the 
molecules,  because  of  the  exceptionally  weak  nature  of  this  bond.  It  seems 
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likely  that  a  multlconfigurat Lon  wavefunction  based  calculation  is  needed  to 
accurately  determine  the  structure  and  atomization  energy  of  the  molecule. 


Table  3.  Results  for  B2 ,  ^2^2 y  anc*  B2H2:  Molecular  Structures, 
Frequencies  and  Electronic  Energy 

Harmonic  SDTQ-MP4  Energy 


Structural 

Parameters 

Frequencies 

(cm-1) 

Atomization  Energy 
(hartrees) 

V 

rR8  I  .910  A 

39  5 

-49  .2456 
0.0469 

B2H: 

rBB  1  *7  iZt  A 

2851 

-49  .9111 

rBH  1*1799 

917 

0.2149 

0BBH  180  *°° 

504 

B2H2  : 

rRB  l  .329  1 A 

2946 

-50  .5843 

rRH  1  .1730  A 

2903 

0 .3906 

® BRH  180 *°° 

678(2) 

685(2) 

We  predict  the  B„H  molecule  to  be  a  linear  free  radical.  We  are  unaware 
of  any  experimental  observation  of  this  species.  The  predicted  BB  bond 
distance  is  much  less  than  that  for  the  B2  molecule,  and  the  predicted 
atomization  energy  is  much  greater. 

The  addition  of  a  hydrogen  to  form  the  ^2^2  mo^ecu^-e  leads  to  a  molecule 
that  appears  to  be  the  borane  analog  of  the  ethylene  molecule,  a  double-bonded 
borane.  The  RB  bond  distance  is  much  shorter  than  that  of  either  the  B2H 
radical  or  the  R?H^  molecule,  and  a  population  analysis  indicates  that  a 
double  bond  connects  the  boron  atoms.  It  seems  likely  that  the  bonding  that 
characterizes  this  molecule  is  a  rarity.  In  addition,  it  should  be  noted  that 
„e  have  found  no  closed-shell,  nonclassical  isomer  of  this  molecule.  This 
distinguishes  the  molecule  from  the  other  Bj  compounds  described  here. 
Schematic  structures  of  these  molecules  are  shown  in  Figure  1  . 
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1  .  Molecular  Structures  for  B2  ,  l^b,  anc* 


1  1 


Figure 


Table  4  summarizes  the  data  for  the  and  B7 molecules.  Both 

classical  and  nonclasslcal  (bridged)  structures  have  been  found  for  each  of 
these  molecules,  and  these  structures  are  shown  in  Figure  2.  Note  that  in 
each  case,  the  '  pron-boron  seperation  is  shorter  for  the  bridged  structure. 
This  consenc'-^e  of  the  three-centered  bonding  in  these  isomers  emphasizes  the 
special  problem  faced  in  developing  a  bond-additivity  correction  scheme  for 
such  compounds.  That  is,  even  though  the  boron-boron  distance  is  shorter  in 
the  bridged  compound,  there  is  not  a  boion-boron  bond.  Electrostatic 
\otential  maps  for  such  molecules  show  that  the  electron  density  is  located 
along  the  three-center  bonds.  The  BAC  should  be  interpreted  in  terms  of  the 
two  three-center  BHB  bonds.  For  both  these  isomer  pairs,  our  results  suggest 
that  the  atomization  energy  for  the  classical  isomer  is  greater  than  that  for 
the  bridged  compound. 

The  structure  of  the  nonclasslcal  Isomers  of  B,Ho  and  B2H4  are 
noteworthy.  All  the  hydrogen  atoms  lie  on  one  side  or  an  imaginary  plane  that 
contains  the  boron-boron  bond.  It  is  interesting  to  conjecture  possible 
transition  states  for  the  interconversion  of  the  two  B2H,(  isomers.  These 
studies,  however,  will  certainly  require  the  application  of  mul t iconf igui ation 
surface  walking  methods. 


Figure  2.  Molecular  Structures  for  the  Classical  and  Nonclassical 

Isomers  of  B^h-j  and 

Data  for  two  isomers  of  and  the  B?H^  molecule  are  presented  in  Table 

S ,  and  the  molecular  structures  are  shown  in  Figure  1.  The  two  isomers 
reported  here  are  clearly  stable  structure's.*  For  the  singly-bridged  isomer, 
the  triangle  formed  by  the  two  borons  and  the  bridge  hydrogen  bisects  the 
angles  formed  by  the  Bl^  groups.  The  doubl v-br idged  has  an  asvmetric  ring 


*It  is  possible  to  suggest  other  bonding  possibilities  for  the  B?H,-  molecule. 
These  include  a  molecule  with  three  bridged  hvdrogens,  or  a  classicallv  bonded 
system.  We  have  recently  determined  the  tri-bridged  isomer  structure. 


Table  4.  Results  for  and  B2H4:  Molecular  Structures, 

Frequencies  and  Electronic  Energies 


Harmonic 

SDTQ-MP4  Energy 

Structural 

Frequencies 

Atomization  Energy 

Parameters 

(  cm" 

-1) 

( har trees) 

B2H3 

Classical : 

rBB 

1  .601  A 

287  1  , 

2748 

-51  .2379 

rBH2 

1  .195  A 

2665, 

1245 

0  .5465 

rBH1 

1.177  A 

1016, 

943 

0  H2BB 

120  .3° 

903, 

450 

G  BBHj 

179  .9° 

434 

Bridged : 

rBB 

l  .487  A 

2937  , 

288  3 

-51  .2370 

rBH 

1  .173  A 

1998  , 

1375 

0  .5456 

r  ! 

BH 

1  .393  A 

1216  , 

841 

G  BRH 

176  .1° 

690, 

634 

°MBH 

1 18  .3° 

38  1 

B2H4 

Cl assical : 

rRB 

1  .686  A 

2729  , 

2725 

-51  .8787 

rBH 

1  .198  A 

2669  , 

2643 

0  .6897 

0  RBH 

121  .6° 

1296  , 

1234 

0  HBH 

1  18  .6° 

1070, 

1070 

869  , 

503 

487  , 

486 

Bridged  : 

rBB 

1  .467  A 

2941  , 

2887 

-51  .8739 

rBH 

1  .172  A 

2144, 

2643 

0  .6849 

r  1 

BH 

1  .365  A 

1415  , 

1383 

G  BBH 

176  .1° 

1257  , 

1  220 

G  HBH 

120.1° 

853, 

503 

G  H,BH 

83  .4° 

566  , 

535 

n 


Table  5.  Results  for  and  I^H^:  Molecular  Structures, 

Frequencies  and  Electronic  Energies 


Structural 

Parameters 


Harmonic 
Frequencies 
( cm- 1 ) 


SDTQ-MP4  Energy 
Atomization  Energy 
( hartrees) 


B2H5  : 


1-Bridge : 


2-Bridge: 


B2H6  : 


rBB 

1  .793 

2848,  2843 

-52  .4526 

rBH 

1  .187 

2718,  2700 

0  .766 

rBH' 

1  .343 

1925,  1339 

0  BBH 

116  .1 

1237,  1195 

0  BH'  B 

83  .7 

1097,  999 

0HBH’ 

126  .7 

933,  853 

726,  497 

478 

rBB' 

1  .77  2 

2815,  2788 

-52  .4451 

rBH' 

1  .347  ,  1  .346 

2707  ,  2265 

0.758 

rBH' 

1  .310 

2073,  1758 

rB'H' 

1  .310 

1677 ,1245 

rBHt 

1  .186 

1073,  983 

rBHc 

1  .188 

966 ,  907 

r B'  H 

l  .184 

851,  832 

°HBH 

122  .4 

408 

Q  BR’  H 

127  .3 

0  B' BHt 

120  .1 

Q 

1 17  .3 

B'  BH 

0  B'  H' B 

83  .7 

rBB 

1  .794 

2823,  2809 

D7.P  Basis 

rRH 

1  .187 

2722,  2702 

-53  .1161 

rBH' 

1  .329 

2238  ,  2010 

0.932 

G  HRH 

1  22  .5 

1875,  1770 

Extended 

0  BH’  B 

84  .9 

1259,  1241 

-53  .1424 

1148,  1063 

0  .9624 

1061,  973 

886,  875 

821,  398 
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that  Is  nearly  planar,  with  terminal  hydrogens  located  in  a  perpendicular 
plane.  The  perturbation  theory  results  predict  that  the  singly-bridged  isomer 
has  a  higher  atomization  energy. 


Figure  3.  Molecular  Structures  for  Two  Bridged  Structures  of 

B7 and  ^2^^ 

We  have  discussed  our  diborane^  results  in  a  previous  paper.  Here  we 
note  that  the  calculations  contain  no  surprises.  We  have  not  attempted  to 
locate  other  possible  stable  isomers  of  this  compound.  Since  B2H^  is  a 
crucial  molecule  In  our  effort  to  devise  a  bond-additivity  scheme  based  upon 
the  results  of  ab  initio  calculations,  we  must  briefly  review  the  results  if 
our  extended  basis  set  calculations  for  this  molecule.  These  results  are 
required,  since  this  molecule  is  the  only  bridged  horane  that  we  have  studi  J 
with  an  extended  basis  set.  The  computations  took  full  advantage  of  the 
molecules  C2v  symmetry.  We  provide  the  extended  basis  set  results  in  Table  . 
along  with  the  OZP  data. 


V.  SUMMARY  OF  B2H6  RF.SULTS 

There  have  been  a  number  of  calculations  dedicated  to  predicting  the 
dimerization  energy  of  horane  and  we  have  reviewed  many  of  those  results  in 
Reference  10.  Our  calculations  employed  the  largest  basis  set,  and  included 
all  the  terms  that  contribute  to  the  electronic  energy  through  fourth-order 
perturbation  theory.  We  cannot  claim  that  the  perturbation  theory  series  has 
converged;  we  rely  on  the  hope  of  judicious  cancellation  of  the  higher  order 
terms  when  we  evaluate  the  atomization  energy.  Additionally,  in  the  case  of 
diborane,  we  have  not  correlated  the  core  orbitals  in  order  that  we  could 
perform  the  fourth-order  calculation  using  the  full  basis  set.  Thus,  in  order 
to  predict  the  dimerization  energy  of  horane,  we  have  also  computed  the 
fourth-order  perturbation  theory  energy  of  BH^  using  the  frozen  core 
approximation.  The  value  we  obtain  is  -43.23  kcal/mole.  This  value  is 


IS 


substantially  lower  than  the  value  reported  by  Ortiz  and  Lipscomb^  in  a 
third-order  perturbation  theory  calculation  (-39.86  kcal/mole)  and  our  value 
of  -40.81  kcal/mole  obtained  using  their  basis  set.  Since  we  know  the 
atomization  energy  of  BH-^ ,  we  can  compute  the  atomization  energy  of  diborane 
directly  using  this  value  of  the  dimerization  energy.  We  obtain  603.9  kcal/ 
mole  (0.9623  h) .  This  is  substantially  lower  than  the  DZP  result,  584.8  kcal/ 
mole  (0  .932  h)  .  We  would,  in  the  spirit  of  the  BAC  approach,  claim  that  the 
energy  difference  is  represented  by  correction  factors  from  four  B-H  bonds  and 
two  B-H-B  three-center  bonds.  We  do  not  consider  a  BB  bond  correction  factor, 
since  there  is  no  BB  bond  in  diborane. 

We  have,  to  this  point,  discussed  our  results  in  terms  of  the  electronic 
energy  differences  alone.  To  estimate  the  more  familiar  heat  of  formation,  we 
must  include  zero-point  energy  differences  and  temperature  effects.  Our 
calculations  predict  a  zero  point  energy  difference  contribution  to  the 
dimerization  energy  of  6.1  kcal/mole,  and  the  temperature  effects  equal  -22.6 
kcal/mole,  for  T=298  K.  Therefore,  we  predict  that  the  dimerization  energy  of 
BH^  is  -39.6  kcal/mole. 

Similarly,  we  can  correct  the  atomization  energy  prediction  for  these 
effects,  and  obtain  the  heat  of  formation  for  diborane,^  1  .6  kcal/mole.  The 
heat  of  formation  for  borane^  equals,  20.6  kcal/mole.  These  estimates  used 
scaled  values  of  the  predicted  harmonic  frequencies  rather  than  experimental 
values.  The  results  described  in  our  previous  paper  used  the  experimental 
frequencies  for  diborane.  While  those  values  are  to  be  preferred  to  the 
scaled  theoretical  data,  similar  data  is  not  available  for  many  of  the  other 
compounds  descrihed  in  this  note.  Consistency  requires  the  choice  of  the 
theoretically  predicted  data. 

VI.  B0NI>- ADDITIVITY  CORRECTIONS  AND  HR  AT  OF  FORMATION  ESTIMATES 

The  data  described  here  is  insufficient  to  accurately  derive  bond- 
additivity  correction  factors  to  provide  the  basis  for  the  development  of  a 
model  thermochemistry  for  the  boranes.  We  have  neither  enough  data  on  the 
boranes  in  general,  nor  enough  accurate  data  upon  which  to  develop  reliable 
BACs .  We  will,  however,  provide  estimates  of  the  heats  of  formation  for  the 
radicals  described  in  this  study  relying  on  the  data  described  for  boranef 3] , 
diborane[6] ,  and  the  results  of  an  extended  basis  set  calculation  for  the 
classical  isomer  of  diborane[4],  where  we  used  the  basis  set  additivitv 

Q 

technique  of  Pople,  et  al .  Using  the  BH^  data,  wc  determine  that  the  BAC  for 

a  boron-hydrogen  bond  is  2.7  kcal/mole,  and  the  diborane  studies  predict  the 

BAC  for  a  B-H-B  three-center  bond  to  be  4.7  kcal/mole,  and  that  of  the  B-B 

single  bond  to  be  7.0  kcal/mole.  We  do  not  have  enough  data  to  estimate  bond- 

length  dependence  od  the  BACs,  nor  do  we  have  enough  data  to  treat  asymmetric 

three-center  bonds  differently  than  we  treat  the  symmetric  bond  of 
£ 

diborane.  The  estimated  heats  of  formation  for  all  the  species  are  collected 
in  Table  6.  These  values  depend  directly  upon  the  assumed  value  of  the  heat 
of  formation  assigned  to  the  boron  atom.  Since  that  parameter  has  a  large 
error  bar,  the  assigned  heats  of  formation  have  a  nontrivial  uncertainty. 

Clearly,  more  work  is  needed  if  one  is  to  develop  a  reliable  method  for 
predicting  the  heats  of  formation  of  the  boranes.  We  must  extend  our  ability 
to  compute  accurate  values  of  the  atomization  energy  to  molecules  containing 
three  or  four  boron  atoms.  Note  that  we  have  chosen  to  describe  the 


correction  factors  in  terms  of  the  bonding  types  that  characterize  the 
boranes,  rather  than  in  terras  of  internuclear  seperations  between  pairs  of 
atoms.  This  selection  is  motivated  by  the  desire  to  develop  group-additivity 
schemes  for  this  class  of  compounds.  With  the  emergence  of  efficient, 
accurate  quantum  chemical  methods  and  the  availability  of  modern 
supercomputers,  one  expects  that  a  sufficient  database  could  he  generated  for 
development  of  an  "empirical"  approach  to  thermochemistry  based  upon  ab  initio 
data . 

Table  6.  Heat  of  Formation  at  298  K  (kcal/mole) 


Molecule 

Hf8 

BH 

103  .5 

BH? 

73  .8 

bh3 

20  .8 

b2h 

176  .3 

»2K2 

1 17  .4 

b2h3(c) 

75  .8 

B2H^(nc) 

73  .3 

b2h4(c) 

38  .4 

B2H4(nc) 

37  .2 

B2H5( lb) 

54  .1 

B2H,-(2b) 

56  .3 

82H6 

1  .6 

VTT.  SUMMARY 

We  have  reported  data  on  the  structure,  harmonic  vibrational  frequencies 
and  heat  of  formation  for  thirteen  mono  and  diboranes.  This  provides  the 
minimal  data  needed  to  estimate  the  heats  of  formation  for  these  compounds 
using  a  consistent  ab  initio  approach.  By  using  extended  basis  fourth-order 
perturbation  theory  methods  to  predict  atomization  energies,  we  develop  the 
basis  for  a  bond-additivity  correction  scheme  that  may  be  used  to 
systematically  correct  calculations  done  with  more  standard  basis  sets.  Since 
the  smaller  calculations  may  be  applied  to  many  more  compounds  than  the 
extended  calculations,  enough  ab  initio  data  may  be  obtained  to  develop  a 
method  for  estimating  thermochemical  properties  based  entirely  upon  ab  initio 
data . 
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